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The direct transition from an insulator to a superconductor (SC) in Fermi systems is a problem 
of long-standing interest, which necessarily goes beyond the standard BCS paradigm of supercon¬ 
ductivity as a Fermi surface instability. We introduce here a simple, translationally-invariant lattice 
fermion model that undergoes a SC-insulator transition (SIT) and elucidate its properties using 
analytical methods and quantum Monte Carlo simulations. We show that there is a fermionic band 
insulator to bosonic insulator crossover in the insulating phase and a BCS-to-BEC crossover in the 
SC. The SIT is always found to be from a bosonic insulator to a BEC-like SC, with an energy gap 
for fermions that remains finite across the SIT. The energy scales that go critical at the SIT are 
the gap to pair excitations in the insulator and the superfluid stiffness in the SC. In addition to 
giving insights into important questions about the SIT in solid state systems, our model should be 
experimentally realizable using ultracold fermions in optical lattices. 


Understanding superconductor-insulator transitions 
has long been an important challenge in condensed mat¬ 
ter physics. The Bose Hubbard model and the Joseph- 
son junction array have led to key insights into the 
SIT in bosonic systems [1-4]. In contrast, there has 
been less progress in understanding the SIT in Fermi 
systems, despite the existence of many electronic sys¬ 
tems that exhibit a direct transition from an insulator 
to a SC. These include, e.g., the disorder-driven SIT in 
thin films [5-7], superconductivity in doped band insu¬ 
lators like SrTiOa [8] and at oxide interfaces [9], and 
the SIT induced in Mott insulators by doping (high-Tc 
cuprates) [10, 11] and pressure (organics) [12]. 

Here we introduce and analyze a two-dimensional (2D) 
translationally-invariant lattice fermion model with local 
interactions that exhibits a direct quantum phase tran¬ 
sition from an insulator to a SC. Our goal is to gain 
insights into a number of key issues in the field of SIT 
through a simple (disorder-free) model that can be ana¬ 
lyzed in great detail, rather than to describe a specific ex¬ 
perimental condensed matter system. We also note that 
our model can be realized experimentally using ultracold 
Fermi atoms in optical lattices. 

In the standard Bardeen-Cooper-Schrieffer (BCS) 
paradigm, superconductivity is a Fermi surface instabil¬ 
ity of a normal metal. The challenge here is to under¬ 
stand how superconductivity arises out of an insulator 
that has no Fermi surface. The key insight from our 
work is that there are Fermi-Bose crossovers in both the 
insulating and in the superconducting phases, and the 
SIT is always between a bosonic insulator and a SC in a 
Bose-Einstein condensation (BEC) regime; see Fig. 1. In 
the weak coupling limit, though, the bosonic regimes on 
either side of the SIT may be narrow. 

Our main results are: 

1) The pairing susceptibility x in the insulator diverges 
and the gap to pair excitations Wpair in the insulator van¬ 


ishes upon approaching the SIT. 

2) The single-particle energy gap Eg for fermions remains 
finite in both phases across the SIT. 

3) The SC state is characterized by a pairing amplitude 
A and a superfluid stiffness Dg, both of which vanish ap¬ 
proaching the SIT. 

4) The insulating state near the SIT is bosonic in the 
sense that the gap to pair excitations is much smaller 
than the fermionic gap. 

5) The SC in the vicinity of the SIT is a BEC-like regime 

with several unusual properties. Its fermionic energy gap 
Eg = [(Eg)^ -I- depends upon both pairing A and 

the insulating gap E^, and its superfluid stiffness Dg is 
much smaller than the energy gap. Its Bogoliubov exci¬ 
tation spectrum exhibits a minimum gap at k = 0 rather 
than on an underlying Fermi surface, giving rise to a gap 
edge singularity that is qualitatively different from the 
BCS result. 

We work with a half-filled attractive Hubbard model on 
a triangular lattice bilayer; the reasons for this particular 
choice of lattice are explained in detail below. Our results 
are based on a variety of analytical approaches, including 
a strong coupling analysis about the atomic limit, a weak- 
coupling analysis of the pairing instability in an insulator 
and mean field theory (MFT). We also present numerical 
results from determinant quantum Monte Carlo (DQMC) 
simulations that are free of the fermion sign-problem. 

Before describing our work in detail, we comment on 
its relationship with the classic paper by Nozieres and 
Pistolesi on “pairing across a semiconducting gap” [13]. 
They used MFT and estimates of phase fluctuations to 
analyze superconductivity in a system with a band gap 
that separates two bands, each with a constant density 
of states (DOS). Building on their ideas, our work goes 
beyond their analysis in terms of what we calculate and 
the methodology used, and this leads to new insights into 
the problem. Our explicit lattice Hamiltonian permits us 
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FIG. 2: Mean-field theory (MFT) results across the T = 0 
SIT in the attractive Hubbard model on a triangular lattice 
bilayer (inset). The SIT is tuned by t/t±_ at fixed filling n = 1 
and attraction \U\/t± = 1.7. The inverse pairing susceptibil¬ 
ity 1/x and gap to pair excitations Wpair in the insulator both 
vanish at the SIT. The SC is characterized by a pairing ampli¬ 
tude A and superfluid stiffness Da, which also vanish at the 
SIT. The single-particle energy gap Eg remains finite across 
the SIT: Eg = Eg in the insulator, Eg = [(Fg)^ + in 

the BEC regime near the SIT and Eg = A is the BCS regime. 


FIG. 1: Schematic phase diagrams at T = 0: (a) in the 
{t, |?7|)-plane at density n = 1, and (b) in the (t, /i)-plane for 
fixed |?7|. The top, middle and bottom insulating lobes in (b) 
correspond to densities of n = 2,1 and 0 (vacuum) respec¬ 
tively. The Fermi insulator (FI)-Bose insulator (BI) crossover 
is defined by the nature of the excitations, single fermion or 
pair, that has the lowest gap. The BCS-BEC crossover in the 
SC is determined by location of the minimum gap in k-space; 
see text. 


to use DQMC and is of a form that can be realized in cold 
atom experiments. We conclude with a comment on the 
connection between our results and other problems - such 
as the disorder-tuned SIT, the superfluid-Mott transition 
for bosons and the BCS-BEC crossover in multi-band 
systems. 

Model: We begin with the constraints on a fermion 
Hamiltonian that realizes a band insulator-SC transition. 
First, we need at least two sites (or orbitals) per unit cell 
to describe a band insulator. Second, we must ensure 
that the attraction needed for SC does not lead to other 
broken symmetries. The attractive Hubbard model on 
a bipartite lattice has an SU(2) symmetry at half-filling, 
with a degeneracy between SC and charge density waves 
(CDW) that leads to Tc = 0 in 2D. To avoid this, we 
choose a non-bipartite lattice. Finally, we want to tune 
the SIT at a fixed commensurate filling. Away from this 
filling the band insulator becomes a metal, and we do not 
get an insulator to SC transition. 


A simple model that meets these criteria is the attrac¬ 
tive Hubbard model on two coupled triangular lattices 
(inset of Fig. 2) with the Hamiltonian 

H = -t H 

(H)||C" {ij)±cr 

- - l/2)(n.4 - 1/2)- (1) 

i<T i 

The spin a =t) i fermion operators at site i are cj^ and 
Cj^, with hopping t between in-plane neighbors and 
t± between interlayer neighbors The chemical po¬ 

tential is /r, the local attraction is |17|, and riia = cl^c^cr- 
Recently the SIT has been studied [14] in an attrac¬ 
tive Hubbard model on a square lattice with near- and 
next-near-neighbor hopping and a staggered (“ionic”) po¬ 
tential to double the unit cell. This model differs from 
ours in that it has one additional parameter and exhibits 
CDW order in a limiting case. More importantly, the 
questions we address and the methodology we use are 
quite different. 

Non-interacting and Atomic Limits: We begin 
with exactly solvable limits in the phase diagram in 
Fig. 1. First, consider the noninteracting [U = 0) sys¬ 
tem with dispersion 


= -t 


\ ' gi(fcx cos sin -^^ ) 


m—1 


t_L cos kz — (2) 
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Here {kx,ky) lie in the triangular lattice Brillouin zone 
and kz takes on values 0 (and tt) for the bonding (and an¬ 
tibonding) band. For fixed n = 1, this implies a transition 
from a band insulator (gap = 2t± — 9t) to a metal 
at t = 2t±/9; see Fig. 1(a). At U = 0 the (t,^) phase 
diagram looks qualitatively similar to Fig. 1(b), except 
that the insulating “lobes” are triangular in shape with 
Mc/t_L = ±1 and the SC is replaced by a metal. (See 
Fig. 6a in Appendix A). 

Next consider the “atomic” limit t = 0, for which the 
lattice breaks up into disconnected (vertical) rungs. We 
solve in Appendix B the two-site Hubbard model on a 
rung for arbitrary t± and \U\. For n = 1 there is a 
crossover from a fermionic insulator to a bosonic insu¬ 
lator at \U\/t_\_~2. For |C/|/tj_ <2, the lowest energy ex¬ 
citation is a single fermion (particle or hole), and hence 
similar to the fermionic band insulator, which is ground 
state for {7 = 0. On the other hand, for \U\/t^_ > 2 the 
lowest energy excitation is a pair of fermions. In the large 
U limit, the ground state is a Mott insulator of bosons, 
with one boson per rung. 

At t = 0 we also find that the extent of the n = 1 
insulating phase, (—/ic,Mc) in Fig. 1(b), is reduced with 
increasing \U\/t±. Thus the n = 0 and 2 lobes grow in 
size relative to n = 1 as \U\/t±_ increases. We note that 
the description of the atomic limit phases as ‘insulators’ 
is justified given that the gap is robust to turning on a 
small tyfO hopping. 

Pairing Instability in the Insulator: We next de¬ 
scribe a weak-coupling theory for the dominant instabil¬ 
ity in the band insulator as we turn on an attraction |{7|. 
The q = 0 dynamical pairing susceptibility in the lad¬ 
der approximation is given by y(w) = [Y7^(a;) — iC/ll , 
with xo(w) = N-'^ “ 2/k)/(2ek - w - f0+). Here 

N is the number of lattice sites, /k is the Fermi func¬ 
tion, and ffk = ^k ~ k-H includes the Hartree shift 
= \U\{n — l)/2. We analyze the problem for n = 1, 
choosing /i in the insulator so that we take a trajectory 
in Fig. 1(a) that goes through the tip of the lobe. (This 
choice of chemical potential is described in Appendix D). 

The divergence of x = a(w = 0) at the SIT is shown 
in Fig. 2. We tune through the SIT by varying 
which controls the gap in the band structure, keeping 
\U\/t± fixed. It is energetically favorable to create pairs 
of particles and of holes when the gain in pair binding 
energy exceeds the band gap. This triggers the SIT, a 
particle-particle channel analog of exciton condensation 
in semiconductors [15, 16]. 

The dynamical pair susceptibility xi^) exhibits a pole 
at the two-particle gap Wpair- We see in Fig. 2 how Wpairj 
the energy to insert a pair into the insulator, goes soft 
and vanishes at the SIT. In the SC, where pairs can be 
inserted into the condensate at no cost, Wpair = 0. We 
show below that the single-particle gap Eg remains finite 
across the SIT. 

Mean Field Theory of SC state: The pairing 


amplitude A = \U\^^{ci^Cil)/N characterizes the SC 
state. We find A and /i from the mean field equations 
l/\U\ = A-iX;ktanh(£;k/2r)/2F;k and (n - l)/2 = 
—^j^(ek/2£'k) tanh(i?k/2T). The Hartree-shifted 
dispersion Sk = £k ~ kH determines the Bogoliubov spec¬ 
trum ifk = + A^, from which we calculate the en¬ 

ergy gap Eg = minAk in the single-particle DOS N{uj). 

The evolution of the T = 0 energy gap Eg across the 
SIT is shown in Fig. 2. In the insulator we call the gap 
Eg = Eg, but once SC sets in, we find that the SC and 
insulating gaps add in quadrature Eg = [{E^Y -f 
For large t/M, the two bands merge, the insulating gap 
Ag collapses and Eg = A, as in BCS theory. 

A single-particle gap Eg that remains finite and a pair- 
gap Wpair that vanishes at the SIT implies that the insu¬ 
lating state close to the SIT is a boson insulator. As dis¬ 
cussed in the Conclusion, these results are very similar to 
those in simple models of the disorder-tuned SIT [17-19]. 

Given that the gap Eg remains finite, what is the crit¬ 
ical energy scale as the SIT is approached from the SC? 
We show in Fig. 2 that the superfluid stiffness Ds goes 
soft at the SIT. Dg is obtained from the q— > 0 limit of 
the transverse current-current correlation function [20]. 

BEC-BCS crossover: We next argue that the SC 
state near the SIT is more akin to the BEG regime than 
to the BCS regime [21]. This is best seen from the ex¬ 
citation spectrum in Fig. 3. For large t/t±, we are in a 
BCS regime, with the usual Bogoliubov dispersion ±ifk 
with weights and v^. The minimum gap A occurs at a 
finite wavevector {kp in weak coupling), which identifies 
an “underlying” Fermi surface (FS) [22]. This gap mini¬ 
mum located on a ID FS contour in 2D k-space leads to 
the well known {E — Eg')~^l'^ singularity in the DOS (in 
addition to van Hove singularities in the band structure). 

The BEC regime near the SIT differs from BCS in a 
variety of ways. The energy gap Eg = [(E^^-l-A^]^/^ is 
located at k = 0. The fact that minAk occurs at a point 
(not a contour) leads to a jumy discontinuity in 2D at the 
gap-edge (not a square-root singularity). This qualitative 
difference in the DOS singularity in the BEC and BCS 
regimes seems not to have been recognized earlier. 

A comparison of the superfluid stiffness Dg and the 
pairing A is also illuminating; see Fig. 2. For large t/t± 
we find Dg ^ A as in BCS theory. Close to the SIT, 
however, we Hnd <C A, a BEC-like regime with well- 
formed pairs but a small phase stiffness. 

Limitations of MFT: The MFT results described 
above give many important insights into the SIT as a 
function of t/t± for fixed \U\/t± < 2, but the approxima¬ 
tions involved give misleading results for large \U\. MFT 
overemphasizes order and predicts a SC state at all t/t± 
for \U\/t± > 2. (See Appendix D for details.) It thus 
fails to describe the Fermi to Bose insulator crossover, 
shown in Fig. 1(a), as well as the SIT that occurs at t/t± 
of order unity for \U\/t± —> oo. 

DQMC results: To investigate the role of quantum 
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FIG. 3: Spectral functions (top) and density of states N{E) (bottom) from T = 0 MFT with \U\/t± — 1.7. The dispersion 
is along the (0,0) to (27r/3, 27r/3) in the Brillouin zone, with line thickness proportional to spectral weight. Left panels: 
BEC regime where the lowest-energy Bogoliubov excitations occur at k = 0, leading to a discontinuity at the gap edge in 
N{E). Right panels: BCS regime where the lowest-energy excitations occur at a finite wavevector ~ Ff, leading to the usual 
inverse-square-root singularity at the gap edge. Dashed lines indicate the Fermi energy. Arrows indicate dispersion minima. 



FIG. 4: DQMC results across the SIT for L x L x 2 bilayer 
systems with \U\/t± = 4 at T = 0.08031^- (a) q = 0 pairing 
structure factor Ps (see text); (b) Superfiuid stiffness Da. 


fluctuations beyond MFT and to obtain quantitative in¬ 
sights at intermediate and large U, we use DQMC simu¬ 
lations [23-25], to solve the triangular bilayer attractive 
Hubbard model Eq. (1). DQMC provides statistically 
exact solutions on lattices of finite size by mapping the 
interacting electron problem onto one of electrons moving 
in an auxiliary field fluctuating in space and imaginary 
time, which is sampled stochastically. (For QMC studies 
of a single triangular lattice, see Ref. 26.) 

We cannot obtain A directly from DQMC, so we com¬ 
pute P, = the q = 0 pairing struc¬ 

ture factor, which equals |A|^ in the infinite-size limit. 
We also compute the superfluid stiffness Dg using the 
standard Kubo formula [20]. We see the SIT in Fig. 4 
from the onset of both |A|^ and Dg as a function of 
t/t± at fixed \U\/t± = 4. We found similar results at 
\U\/tj_ = 3 (not shown). 

We see that Dg increases monotonically with t/t± but 
Ps exhibits non-monotonic behavior, similar to the MFT 
results for A, which we can understand as follows. In 
the BEC regime, close to the SIT, the order parameter 
A increases with t/t±, however, eventually the increase 
in bandwidth leads to a smaller normal-state DOS, and 
the BCS A decreases. 

Finally, we show in Fig. 5 that the persistence of a finite 
single-particle gap Eg from the insulator to the SC is not 
an artifact of MFT, and is clearly seen in the DQMC 
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FIG. 5: DQMC density of states N{uj), calculated using 
the maximum entropy method, clearly shows the persistence 
of the single-particle energy gap across the SIT from the insu¬ 
lator to the SC. The dashed black line indicates the SIT. (a) 
False color plot of N{oj) as a function of t/t± on a 12 x 12 x 2 
bilayer with \U\/t±^ = 4 and T = 0.0803t^; (b) N(ijj) for 
specific values of t/t±. 


results for the DOS. The DOS N{lo) was obtained from 
analytic continuation of DQMC data using the maximum 
entropy method [27]. A detailed analysis of our DQMC 
results will be presented elsewhere. 

Conclusions: Our main results are summarized in 
the Introduction. We conclude by discussing the parallels 
between our results and several problems of current inter¬ 
est. First, we comment on the very different problem of 
the disorder-driven SIT in the attractive Hubbard model, 
analyzed using a spatially inhomogeneous Bogoliubov- 
de Gennes MFT [17, 18] and DQMC [19, 28, 29]. The 
predicted [17, 18] persistence of the single-particle gap 
across the SIT has been seen in STM experiments [30- 
32]. The mechanism leading to the insulating gap cru¬ 
cially involves inhomogeneity and an unusual insulator 
with localized pairs, and is quite different from the gap 
that persists across the SIT analyzed here. 

The collapse of the two-particle gap Wpair on the in¬ 
sulating side of the SIT and that of the stiffness Dg on 
the SC side also have interesting parallels between the 
disorder-driven SIT [19] and the SIT studied here. How¬ 
ever, the disorder-driven SIT is hard to analyze analyti¬ 


cally in the manner presented here given the complexity 
of treating both disorder and interactions. 

Next, we note some similarities of our results which 
are on a fermionic system with the superfluid-Mott tran¬ 
sition in the repulsive Bose-Hubbard model (BHM) [1]. 
Of course, the critical behavior at the SIT in our model 
is expected to be in the same universality class as that 
in the BHM. What is much more interesting, however, is 
that the phase diagrams of the two models in the (/i, t)- 
plane are so similar - see Fig. 1(b) - even though, for 
small |C/|, our fermionic system cannot be mapped onto 
bosons. 

Finally, we comment that the BCS-BEC crossover 
problem [21], investigated in great detail for ultracold 
Fermi gases, has almost exclusively focussed on the case 
of a single band. The BCS-BEC crossover described here 
has interesting new features that relate to the two bands 
and the underlying insulating band gap. The study of 
the BCS-BEC crossover in multi-band systems is in its 
infancy and there are indications that these ideas may be 
relevant to SC in materials such as FeSeTe [33]. 

Acknowledgments: We gratefully acknowledge sup¬ 
port from NSF DMR-1410364 (MR), DOE DE-FG02- 
07ER46423 (NT), and from the UC Office of the Presi¬ 
dent (CC, RTS). 


Appendix A: Non-interacting limit 


First consider the non-interacting tight-binding model 
on a triangular bilayer (Eq. (1) with \U\ = 0), whose 
dispersion relation is given by Eq. (2). The density of 
states (DOS) is 


9iE) = - 


5tri 


E — t± 


9tr 


E -\- t_L 
t 


(A3) 


where gtriis) = — Im Gtri(e+t0+) is the density of states 
of a single triangular lattice, and Gtri is the triangular 
lattice Green function [34], 


GtrK^) — 


dpdq 


1 


0 TT^ £ — 2(cosp-f cosq-I-cospcosg) 
1 



2 (1-1) (1-1)'^'(I+3)'^^ 2^ 


(A4) 


where K is the complete elliptic integral of the first kind 
as implemented by the Mathematica function EllipticK. 

For small t, this model possesses a valence band with 
energies —t± — 6t < E < —t± + 3t and a conduction 
band with energies -|-t_L — 6t < E < -|-t_L + 3t. When 
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t > 2tj_/9 the bands overlap to form a single band. The 
model behaves as a band insulator or metal depending on 
where the Fermi energy /i lies within the bandstructure. 
This leads to the phase diagram in Fig. 6a. 

Appendix B: Atomic Limit 

Now consider the triangular bilayer Hubbard attractive 
model in the “atomic limit” {t = 0, \U\ > 0). The system 
now consists of independent two-site Hubbard models on 
disconnected rungs of the bilayer, each rung described by 

H = -t±Y^ (cVcfia + ^-c-) 

<7 

- \u\ - l/2){nii “ 1/2) (A5) 

icr i 

where i = A,B distinguishes the two sites. Each two-site 
system can be occupied by iV = 0,1,2,3, or 4 fermions. 
Exact diagonalization in the basis of 16 Fock states shows 
that the lowest-energy state in each subspace of fermion 
number N is 

Eo = -^\U\+2^i 
El = — 1± + /i 

E2 = -^l\Uf+4tl 

Es = -t± - ^ 

E^ = -^\U\-2^i. (A6) 

Fig. 7 shows the energy levels as a function of U (for 
/X = 0). The ground state is always at half-filling {N = 2, 
black curve). 

For \U\ /t± < 2, the lowest energy excitation out of 
the N = 2 subspace is a fermionic excitation, i.e., adding 
or removing one fermion to get to iV = 1 or = 3 (red 
line). In this regime the model might be said to be a 
“Fermi insulator,” with a single-particle gap 

Eg = min (Es - E 2 , Ei - E 2 ) ■ (A7) 

For \U\ /t± > 2, the lowest energy excitation out of the 
N = 2 subspace is a bosonic excitation, i.e., adding or 
removing a pair of fermions to get to the N = 0 ov N = 4 
subspace (blue line). In this regime the model is a “Bose 
insulator,” with a two-particle gap 

Wpair = min {Ei - E 2 ,Eq- E 2 ) ■ (A8) 


Appendix C: Atomic Limit plus Hopping 

We now consider the effects of a small in-plane hopping 
t ^ t± on the atomic limit results. 


Corrections to ground state: For t = 0, the half- 
filled ground state of the system can be schematically 
written as |... 2222 ...). Here each “2” represents the 
fact that there are two fermions on every rung in the bi¬ 
layer lattice when the density per site is n = 1. In fact, for 
each rung, ” 2” is a specific linear combination of six states 
Iti, 0), |0, ti), |ct, a') with two fermions on one rung, with 
definite amplitudes that depend on \U\/t±. In the pres¬ 
ence of finite in-plane hopping t > 0, the ground state 
develops an admixture of states such as |... 2312 ...) and 
|... 2402 ...). We will ignore these corrections. 

Corrections to excited state with one extra par¬ 
ticle: If t = 0, the lowest single-particle excited state is 
of the form |... 2322 ...), and the single-particle gap is 

Eg{0) =E3-E2 = \/l\Uf + Atl -t±-E (A9) 

If t > 0, the lowest single-particle excited state is 
mainly a superposition of states such as |... 2322 ...) and 
|... 2232 ...), connected by hopping t. Thus we might 
expect that the new single-particle gap is 

Eg{U)^ Eg-at (AIO) 

where a is a constant of order unity. 

In the strong-coupling limit \U\ ^ tj_. Eg « |C7| /2, so 
we expect the single-particle gap to collapse to zero at 
about t ~ |17|. However, this is not an insulator-metal 
transition (as we see in the next section). Rather, the 
collapse of the “nominal” single-particle gap represents 
a crossover from a BEC superfluid to a BCS superfluid 
regime. This behavior, Ibec-bcs oc \U\, is shown as a 
dashed line in Fig. 1(a) of the text. 

In the weak-coupling limit |C/1 ^ t±, the above approx¬ 
imation suggests a transition at t ~ tj_. This is a crude 
estimate; we actually know that the transition actually 
occurs at t « 2t±l9 as \U\ vanishes, as shown earlier. 

Corrections to excited state with two extra par¬ 
ticles: If t = 0, the lowest two-particle excited state is 
of the form |... 2422 ...), and the two-particle gap is 

Wpair(O) =Ei-E2 = i |C/| - 2/x. 

(All) 

If t > 0, a local two-particle excitation can hop from one 
rung to an adjacent rung in a two-step process. The ef¬ 
fective boson hopping is of the order of tboson = 4f^/ |[/|. 
Hence the two-particle gap is reduced to 

^pair(If) ~ ^pair(O) O tb oson (A12) 

where a' is a constant of order unity. In the limit |C7| ^ 
t±, we have Wpair(O) « 4tj_^/|17|. This suggests that 
the two-particle gap Wpair falls to zero when 4t±^/ \ U\ — 
a X 4t^/|{7| = 0, i.e., when t ~ t±. Then the system 
becomes overrun by two-particle excitations and makes 
a transition from BI to BEC. This estimate is crude, as it 
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(a) (b) (c) 


FIG. 6: (a) Phase diagram of the non-interacting triangular bilayer tight-binding model. As t increases, the valence band 

and conduction band broaden and eventually overlap, causing a n = 1 band insulator to metal phase transition, (b) Phase 
diagram of the attractive Hubbard triangular bilayer as function of /iefi at |f/|/fx = 1.7, where includes the Hartree shift. 
The triangular regions shrink to lobes at the locus of the pairing instability. A portion of the insulating phase in (a) phase is 
converted into a BEC superfluid, whereas the metallic phase turns into a BCS superfluid, (c) Phase diagram of the attractive 
Hubbard triangular bilayer as function of actual chemical potential /r at \U\/t± = 1.7. The n = 2 insulating lobe is displaced 
downward by \U\ /2, whereas the n = 0 lobe moves up by \U\ /2. In the superfluid, lines are distorted due to Hartree shifts by 

non-integer n. 



FIG. 7: Energy levels for a two-site Hubbard model as a 
function of U, for fj, = 0. 

neglects a large number of corrections to the initial and 
final state. Nevertheless, we expect that the general idea 
is correct, i.e., for \U\ 3> t± the BI-BEC phase boundary 
in the (t, U) plane has a vertical asymptote at some finite 
value of t/t±. This behavior is illustrated in Fig. 1(a) of 
the text. 

Appendix D: Pairing Susceptibility and Mean-Field 
Theory 

We next present a mean-held approach, which works 
best in the \U\ <C t regime, and thus complements the 



f/fj. 

FIG. 8: Schematic phase diagram from mean-field theory. 

atomic limit analysis presented in the previous two ap¬ 
pendices. 

Pairing instability: As discussed in the text, we hrst 
calculate the bare particle-particle channel susceptibility, 
with center-of-mass momentum q = 0, given by 

k 

Here Ek = ~ takes into account the Hartree shift 

= \U\{n— l)/2 with the bare dispersion given by 
Eq. (2). n is the density, N the number of lattice sites 
and /k the Fermi function. The k-sum is over the two 
bands /cj, = 0, tt and over the 2D Brillouin zone {kx, ky). 




















In the metallic phase xo has the well-known Inw di¬ 
vergence of BCS theory. However, we concentrate here 
on the pairing instability in the insulating phase, where 
Xo is finite. The T = 0 result for the n = 1 insula¬ 
tor, which corresponds to a choice of chemical potential 
—t± + 3t < fj, < t± — 6t, can be shown to be 


Xo(w) 


1 

4t 


G, 


tri 


^w/2-|-/i-|-t_L^ 


^ w/2 -|- /X — t_L ^ 

(AM) 


We are able to write the susceptibility xo in terms of a 
single-particle Green’s function because we are looking 
only at pairs with total momentum q = 0, built up from 
two single-particle excitations with momenta ±k, and 
= ^-k- hnd that similar expressions hold for 

the n = 0 and n = 2 insulating states. 

We then calculate the pairing susceptibility in the lad¬ 
der approximation, 


x(w) 


1 

Xo'M-lGl' 


(A15) 


We find the critical in-plane hopping tc in MFT by lo¬ 
cating the value of t at which xo(w = 0) = 1/|C/|, so that 
x(a;) diverges. 

Figure 6c shows the MFT phase diagram in the 
plane, for a fixed coupling \U\ ltj_ = 1.8. In the insulating 
phase there is a finite pairing susceptibility xoj which can 
lead to a pairing instability in the presence of attraction 
\U\. As \U\ increases, more and more of the insulating 
lobe is converted into a BEG superfluid. 

Figure 8 shows the MFT phase diagram in the {t, U) 
plane at half-filling. As |I7| increases, the system under¬ 
goes a transition from an insulator (INS) to a superfluid 
(BEG), due to the pairing instability. At the end of this 
Appendix we will discuss the limitations of MFT and 
why the MFT results in Fig. 8 look so different from the 
more correct results shown in Fig. 1(a). 

Ordered state: In the superconducting state, there 
is a finite pairing amplitude A = {\U\/N) X]i(ciT® 4 )- We 
solve the mean-field theory (MFT) equations 

1/|{7| = A-i^tanh(Ak/2r)/2£;k, (A16) 

k 

(n - l)/2 = -N-^ ^(ek/2Ak) tanh(£;k/2r). (A17) 

k 


to find A and /x, where the Hartree-shifted dispersion 
£k = £k ~ determines the Bogoliubov quasiparticle 
spectrum Ak = \Js\^ -I- A^. Using this we can calcu¬ 
late the energy gap Eg = minEk in the DOS N{uj) for 
single-particle excitations. (The k-integrals over the 2D 
Brillouin of the triangular lattice are performed by reduc¬ 
ing them to ID integrals using the exact density of states 
(DOS) in terms of the triangular lattice Green function.) 
The BEC-BGS boundary (dotted black curves in Figs. 6b 


and 6c) is determined by the criterion that the disper¬ 
sion minimum occurs at k = 0 in the BEG regime and 
at k 7 ^ 0 in the BGS regime. 

Choice of chemical potential: We are interested in 
studying the superconductor-insulator transition tuned 
by hopping t, rather than by fermion density n. Thus, 
in the superconducting phase we choose the chemical po¬ 
tential ^ (according to Eq. (A17)) such that the average 
density corresponds to half-filling, n = 1. In the insu¬ 
lating phase, we choose /x such that the two-particle ex¬ 
citation gap is particle-hole-symmetric. In other words, 
the energy cost of adding a pair of fermions, Wea, is equal 
to the energy cost of removing a pair of fermions, ujh 2 - 
This choice of /x corresponds to the dashed blue curve in 
Figs. 6b and 6c, which bisects the n = 1 insulator lobe, 
and passes through the tip of the lobe. The quantities in 
Fig. 2 of the text are plotted for this choice of 

Note that the single-particle excitation gap. Eg, is not 
particle-hole-symmetric along this trajectory. The value 
of Eg plotted in Fig. 2 of the text corresponds to the 
smaller of the particle gap and the hole gap, min(ife. Eh). 

If one attempts to choose /x to make E^ = Eh, one 
finds that uje 2 7^ w/jj. This /i(t) trajectory exits through 
the side of the lobe rather than the tip of the lobe, repre¬ 
senting a density-tuned transition rather than a hopping- 
tuned one. 

Strengths and Limitations of MFT: MFT allows 
us to track the behavior of many quantities (see Fig. 2 in 
the text) that are intuitively meaningful, but not accessi¬ 
ble in DQMG or other methods. For fixed U and varying 
t, MFT gives a good idea of the general behavior of these 
quantities. 

In MFT the interaction is decoupled exclusively in 
terms of the order parameter. Hence MFT fails to cap¬ 
ture important correlations in the insulating state. In 
the atomic limit {t = 0), as |I7| increases, the ground 
state wavefunction adjusts itself to take advantage of the 
attraction to produce a binding energy. There is no su¬ 
perfluid state at t = 0 (as would have been predicted by 
MFT). Gompare Fig. 8 and Fig. 1(a). MFT gives bad 
results as a function of U. 

In MFT, the insulator-superfluid boundary \U\ (t) is 
a decreasing function, whereas atomic limit and DQMG 
results suggest that it is an increasing function as shown 
in Fig. 1(a). Thus we can no longer think of the insulator- 
SC transition as due to a pairing instability triggered by 
increasing \U\. Nevertheless, we can still understand the 
insulator-SC transition at fixed U, where a change in 
bandstructure (e.g., reduction in band gap) that causes 
Xo to increase beyond 1/|C/| and precipitates a pairing 
instability. 
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